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We investigate screening effects of the medium on the potential interaction between two static 
'charges' for different models of dense plasmas in the one-boson exchange approximation. The 
potential can exhibit an oscillatory behavior, which is related to the analytic structure of the corre- 
sponding boson propagators in the complex q-plane. We have first revisited the one-pion exchange 
in a nuclear medium. In addition to Friedel oscillations, which are associated to branch cuts in the 
", g-plane, there appears another oscillatory component, which arises from a pole on the pion propa- 

■ gator. This pole is located appart from the axes, giving rise to an oscillating Yukawa-like potential. 

' Therefore, we call this phenomenon 'Yukawa oscillations'. This phenomenon does not appear in 

the Debye component of the QED screened potential, even if the coupling constant is artificially 
> | ■ increased. We have also studied a model of QCD quark-gluon plasma. In this case, the one-gluon 

propagator also shows this kind of poles. At high densities and/or temperatures, where one expects 
perturbative QCD to be valid, the pole shifts towards large momenta. 
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I. INTRODUCTION 



' The natural framework for the analysis of high-density matter under conditions where relativistic effects 

are essential is Quantum Field Theory Q, where dense plasmas are described (at fundamental or phe- 
nomenological level) as assemblies of fermions (electrons, nucleons, quarks, etc..) interacting through 
the exchange of bosons (photons, mesons, gluons, etc.). The dynamics of these systems is defined by 
Lagrangian models whose solution (in general, in a given approximation), allows for the calculation of 
macroscopic properties. One interesting issue which can be obtained from the Lagrangian (using, for 
example, the one-boson exchange approximation " OBE" ) is the description of the two-body interaction 
in vacuum, in terms of a phenomenological potential. Simple instances are the Coulomb and Yukawa 
potentials associated, respectively, to the exchange of massless and massive bosons. 

A consequence of the collective behavior in a plasma is the appearance of screening effects on the two-body 
interaction. In a classical electromagnetic plasma, the Coulomb-like interaction between two static charges 
in vacuum becomes, when modified by the medium, the exponentially damped Debye interaction. How- 
ever, in degenerate quantum electromagnetic plasmas the in-medium potential becomes long-ranged and 
oscillatory, both in the non-relativistic g and relativistic H, j| cases. These so-called Friedel oscillations 
have observable consequences ||. 

The phenomenological description of the nuclear interaction in vacuum through the exchange of massive 
mesons Q leads, in the static OBE approximation, to an interparticle potential which is exponentially 
damped with the distance @ . The effects of the medium in a degenerate nuclear plasma, as obtained in 
the same approximation Q, or in a direct calculation of the interaction energy ||, introduce qualitative 
modifications on the two-particle interaction potential, which becomes long-ranged and oscillatory, as in 
the case of the degenerate electromagnetic plasma. The analysis of the screened potential in both cases 
allows for a separation into two components: a Friedel-like component, which is damped as a power of 
the distance and dominates the long range, and an exponentially damped Yukawa-like (or Debye-like in 
QED) component. We must stress, however, that there are qualitative differences between this Yukawa 
component and the Debye screening of QED: the first one is, in fact, oscillatory below some critical 
temperature, whereas not such oscillations are present in the Debye component in the QED plasma. 
Friedel oscillations, which are present in both cases, fade away for lower temperatures ||. As we will 
discuss later, these Yukawa oscillations are also present in other physical systems such as a quark-gluon 
plasma, for which Kapusta et al. M] have been evidenced the existence of Friedel oscillations. 



: In the standard terminology, this approach is known as 'Quantum Hadrodynamics' or QHD Q. 
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The analytic structure of the dressed propagators governs the behavior of the screened potentia l [B |. From 
the mathematical point of view, Friedel oscillations are associated to the Kohn singularity in the 
matter part of the zero-temperature polarization tensor, which shows branch cuts starting at a transferred 
momentum q = ±2p/, where pf is the Fermi momentum of fermions in the plasma. On the other hand, 
Yukawa oscillations are associated to a pole of the propagator in the complex q -plane which has a non- 
vanishing real part, in contrast to Dcbye screening, which arises from a pole on the imaginary axis. 

In this paper we analyze the two-particle interaction in dense plasmas, using the above mentioned models. 
Here we are not interested in a precise quantitative description of actual dense matter using more or less 
adapted theoretical models. Our main concern is rather to evidence the above mentioned characteristics 
of the screening, as a result from basic properties which seem to be present in some relativistic models 
describing dense matter and, eventually, to discuss the limitations of such models. 

The ground states of the plasma in the electromagnetic and in the quark-gluon cases are described as 
relativistic free Fermi gases in thermodynamical equilibrium. For the QHD nuclear plasma, the ground 
state is described by the relativistic Hartree approximation (RHA) jl2| , which defines the thermodynamics 
at finite temperature and includes the contributions of vacuum fluctuations. The propagators of the 
exchanged bosons which are necessary to obtain the screened potential are dressed by the medium. Their 
expressions have been obtained in a self-consistent one-loop approximation, where matter polarization 
contributions were calculated at finite temperature and vacuum polarization contributions were included. 
Such approximation, for QHD, is the first order of a kinetic cluster expansion around the relativistic 
Hartree ground state Jl^|l4| and coincide at T = with the usual one- loop approximation jl5| . 

This paper is organized as follows. In section 2 we introduce the essential equations which are necessary 
for the analysis of the screening in the models. In section 3 we study the analytic structure of the meson 
propagators and the oscillatory behavior of the screened potentials. We conclude in section 4 with a 
discussion on the eventual manifestations in actual systems of the phenomena found here. 
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II. BASIC EQUATIONS 



The static two-particle interaction potential in the medium can be obtained from a linear response analysis, 
and can be put under the general form (!(]] : 

vl , Q1Q2 . r , qe l * r 

47T 2 r J_ 00 H q 2 + F(uj = Q,q) 



Here, r is the distance between the positions of the interacting particles, Q\ and Q2 their corresponding 
"charges" and k = (u),q) the transferred four- momentum. A similar expression can be derived from the 
OBE approximation in the static (co = ) limit ||. Then, G(k) — _ k 'i^ F ^ corresponds to the boson 
propagator, and F(k) is related to the boson self-energy. 



To be more explicit, we need to consider each case separately. As representative of QHD, we will study 
the one-pion exchange. In this case, one has to replace Q1Q2 — > — g 2 , with g^ the pion-nucleon coupling 
constant, and take into account the pion mass. This means that the pion propagator G^(uj = 0, q) must 
be written as : 



G„(lj = 0,q) = - 



1 



q 2 + ml - g^m a < a > -g%U n (q) 



where II 7r (fc) is the pion polarization, which can be decomposed into a matter and a renormalized vacuum 
contribution U n (k) = II™ at (fc) + II™ c (fc) and m w is the pion mass. We have included a pion-sigma 
coupling term in the Lagrangian, as in ||, with g aTi = 3.77. Finally, m a = 550 MeV represents the a- 
meson mass and < a > the self-consistent mean-field value of this meson field. The renormalized vacuum 
contribution to the pion polarization is given by: 



rC c (fc) = (2^)- 2 {- a fc 2 iog 



1 + a 
l-o 



m 2 + 2 (to 2 - M 2 ) 



M 2 1 

2(fc 2 - 2Af 2 )log(— )- 4(m - M) 2 (l + an 2 ) - (— )arctan(— ) 
to ao ao 



2(to - M) 2 (l + a 2 ) + fc 2 (l 



2m 2 



2M< 



where a = y/l — AM 2 /k 2 , ao = -\/4to 2 /to 2 — 1 and M is the nucleon effective mass, calculated in 
the Hartree approximation. In the next section we give the expression for the matter contribution. 

From the expression (|l]) we see that the analytic structure of the propagators on the u> — axis becomes 
the essential magnitude which governs the behavior of the screening. In the QHD case there are poles of the 
propagator on this axis associated to the "tachyonic" branches, which come from the vacuum polarization 
terms [ piS] , Such poles arise at large momentum transfer, where the internal structure of the nucleon 
invalidates the point-particle approach of the model and thus, the effect of this pole on the calculated 
in-medium potential is spurious. We can try to account phenomenologically for the nucleon structure 
by the introduction of monopolar form factors ff(k) on the vertex of the loop expansion and the OBE 
diagrams. To do this we have to replace, in the previous equations 

gl - glff(k) 
ff(k) = (A 2 - ^ 2 )/(A 2 - fc 2 ) 



where is a cut-off parameter. When the value of this parameter is properly chosen, the spurious "tachy- 
onic" branches (and the associated pole of the propagator on the to = axis) disappear. Of course, one 
could also eliminate the poles associated to these branches by simply discarding the renormalized vacuum 



polarization contribution. However, this contribution is essential in obtaining the physical propagation 
modes. If neglected, new unphysical space-like modes (which introduce new poles in the propagators at 
small q) are present and can not be eliminated by form factors The standard expression of the 

form factors Eq. (^) fits well the known data in vacuum [M. At finite density and temperature the form 
factor is expected to be modified by the medium. In a simple bag- model picture the nucleons "swell" as 
density increases (EMC effect), due to a reduction of the bag pressure. Higher-order calculations of the 
vertex function p7[ , where the surrounding meson cloud defines the size of the nucleon, leads to a similar 
behavior. Here we shall keep the form factor unchanged in the medium according to our level of approxi- 
mation, but it must be emphasized that, in an improved calculation, the above mentioned effect will cause 
a stronger suppression of high-momentum exchanges at larger densities. Nevertheless, these modifications 
would affect the short-range part of the interaction, and will not modify our main conclusions, which 
concern the long-range behavior of the screened potential. 

We will also analyze the one-gluon exchange screened potential on a color-singlet state of a quark-antiquark 
pair. For this analysis, we have considered a SU{N) quark-gluon plasma with Nf quark flavors. One has 
to replace 

n n (N 2 - 1).9 2 r ,, 
Q1Q2 -> ( ' 

in this case in Eq. ([l]), where g will be the coupling constant. In order to simplify our analysis , we 
have restricted ourselves to a plasma with only u and d massless quarks with the same chemical potential. 
The denominator in Eq. (Q) is then related to the one-loop gluon polarization in the static limit. The 
corresponding formulae can be found, for example, in fl8ll , fl9| and [ppf . We will write explicit expressions 
for zero temperature in the next section. 



III. ANALYSIS OF THE IN-MEDIUM INTERACTIONS. 



As discussed in the previous section, in the case of the pion- nucleon plasma, one has to evaluate Eq. (Q) 
with the quasi-pion propagator given by Eq.(^J) and the introduction of the form factors given by the 
replacement (|4|). The integration can be transformed to a contour integral in the complex g-plane by 
analytical continuation of H(q) = H(ui — 0,q). Let us first discuss the situation at zero temperature. We 
will summarize the main results found in [01. The matter contribution can be written as: 



nr 4 (<z) = -w* 

+aq 2 log 



4s fPf - (4M 

( a Pf + £ f 
\apf-ef 



2q 2 ) log 



M 



2e f q log 



2p f 



2 Pf 



(6) 



Now a = v/l + 4M 2 /q 2 (we have corrected two misprints appearing in Eq.(4.8) of reference |9(]). One 
has on the upper-half plane the analytical structure of the pion propagator shown in Fig. la)0, where we 
have also shown the contour path to be used in evaluating the integral Eq. . The branch cuts starting 
at q = ±2pf appear from the logarithms in Eq.(^) and have their origin on the sharp form of the Fermi 
surface in momentum space at T = 0. They give rise to Friedel oscillations. For the nuclear plasma, they 
could manifest as an instability of the uniform configuration in the Hartree approximation, thus allowing 
for the possibility of new phases in nuclear matter J22|. Alternatively to the use of contour integrals, 
one can integrate by parts Eq. (|l|) , following a systematic method given by Lighthill |2j| ; then one finds 
that H™ at (q) has a pole on the second derivative : the Kohn singularity In addition to this, one has 
the pole contributions at the points q — ±q r + iqi shown in Fig. la) as the small crosses with a circular 
contour around them, which we will refer to as the Yukawa pole. Then the potential V(r) , for large r will 
consist on two terms : 



2 Aside the branch cuts at q — ±2p/ , there are also branch cuts at q = ±2iAf (not shown in the figure). They correspond to 
the thereshold of particle-antiparticle production. However, they will contribute only for the short-distance range (r < 1/m). 
Here, we will neglect this contribution. 
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V(r -f oo) ~ V^(r) + V Y {r) 



The first component Vf{t) is the Fricdel component, which is oscillatory and decays as a power-law of r 
for large distances, while Vy(r) is the oscillatory exponentially-damped term (Yukawa term), which has 
the form : 

V Y (r) = A- sin(q r r)e~ qir 



where the amplitude A, the "effective" pion mass qi and the "frequency" q r depend on the plasma density. 
At low densities, the frequency q r becomes imaginary and the oscillation disappears. When density vanishes 
this component reduces to the usual Yukawa potential 0. 



We will now discuss the situation for finite temperature. At non-zero temperature T, the Fermi surface 
is spread over a thickness T and the Kohn singularity disappears. However, at small temperatures the 
second derivative in the polarization T) keeps still an important jump around q = ±2p/. This 

can be seen in Fig. 2, where we have plotted the second derivative of the polarization as a function of 
q/m, for various temperatures at saturation density (which corresponds to pf ~ 0.3m ). We see the Kohn 
singularity at T = and the evolution of the jump around q — 2pf as temperature increases. For 
temperatures around 20 MeV the propagator is a smooth function of q. 

At low but finite temperatures ((e/ — M)/T 3> 1), with ej the nucleon Fermi energy, the Friedel potential 
becomes 

V F ir,T) w B{Tp 2 f /r 2 )/v 2 (2p s ,T) cos(2 P/ r)/ sinh( 2lTrT£ f ) 

p f 

where 

u(q,T) = (q 2 + Al)(q 2 + ml) - g*(A£ - ml)U T (q,T) 

and B is a constant. In the limit T = , Eq.(^) behaves as 

V F (r) » cos(2P f r)/r 3 



As temperature increases, the Friedel component is exponentially reduced with T (see Eq.(Q)). In the 
present pion-nucleon model this component fades away, at saturation density, as temperature approaches 
T ~ 20MeV. 

The oscillations associated to the Yukawa component are present up to T ~ 60MeV. We have shown this 
behavior in the next figures, where we plot the real part q r (Fig. 3) and the imaginary part (Fig. 4) 
of the Yukawa pole, both in units of the free nucleon mass m as a function of the temperature, for three 
values of the density, defined by a given pf at T = : Pf = 0.3, corresponding to saturation density (solid 
line), pf = 0.35 (dotted line) and pf — 0.4 (dashed line). As it is readily seen, the pole shifts towards the 
imaginary axis as temperature increases, and remains as a pure imaginary pole for temperatures above 
some critical value T c , which depends on the density. As a consequence, Yukawa oscillations are present 
only if T < T c . This transition can also be presented, from a formal point of view, as q r 'becoming 
imaginary' (and negative), thus being effectively substracted from the value of qi and producing the kink 
observed on Fig. 4. 

Now the question is wether the Yukawa pole found in the pion propagator is a peculiarity of this model 
or not. We have performed a similar study for the screened photon propagator in QED and the screened 
gluon propagator in QCD. In the first case, we found a pole on the imaginary axis for all the values of 
the density and temperature we considered. Of course, this imaginary pole corresponds to the well-known 
Debye screening phenomenon. An important difference with QHD is the strength of the coupling constant 

: in QED one has a = y^, whereas for the pion-nucleon model it is much larger (fj- ~ 14). We have 
artificially increased the value of the QED coupling constant, by taking values a ~ 10 in order to test if 
the Yukawa pole arises in QHD as a consequence of the large value of the coupling constant. However, 
even in this case, we have not found such pole in the electromagnetic case. Therefore we conclude that 



the strength of the coupling constant is not the main reason for the appearance of Yukawa oscillations, 
and they must be instead related to the nature of the model. 



The situation is completely different in QCD. Let us consider the gluon propagator at the one- loop ap- 
proximation, for the model described in the previous section. The denominator in Eq.(|l|) can be written 
as : 

D{q) = q 2 + F(q) = q 2 + F ma \q) + F vac {q) 



where F mat (q) and F vac (q) stand for the matter and vacuum contributions, respectively. We follow the 
formulae given in Jicfl to write these contributions (calculations in this reference are made using the 
temporal axial gauge). At zero temperature, the matter term is given by 



F mat (q) 



48tt 



M + 2 f (V, - 3, 2 ) log (|i±i) - 2q 2 log ( 



q 2 — Ap 2 f 



and the remaining terms can be combined to give 



f + F vac {q) = JL-(11JV - 2N f )q 2 log ( ^ 



QCD 



Using Eqs. (|l|) , (PH) , (|l^) and factoring out by defining 



D(q) 



48tt 



one finds: 



V(r) = - "K^^m r 



Nr 



/(?) 



Thus, the potential V(r) is independent of g 2 . At zero temperature, the analytic structure of the gluon 
propagator is shown in Fig. lb). Aside the branch cuts at q = ±2pf , there is also a branch cut on the 
imaginary axis. 

At zero density (p/ =0), F mat (q) vanishes, and one has 



f(q) = (UN~2N f )q 2 \og 




Therefore, the gluon propagator has poles on the real axis (at q — ±Aqcd )j an d the long-distance 
behavior of V(r) is dominated by these poles. We have analyzed the behavior of this pole as pf increases, 
for a chosen value of Aqcd = 200 MeV. First, the position of this pole shifts from the value of Aqcd, and 
a second pole at smaller q appears. However, for pf > 120 MeV both of them disappear and 1/D(q) will 
not have poles on the real axis. Since the deconfinement transition is thought to take place (at T — 0) at 
a density such that pf > 300 MeV, we can assume that no poles on the real axis will be present in the 
quark-gluon plasma phase. Instead, a pole with both a real and an imaginary part appears at qo = q r + iqi, 
as shown in Fig. lb). This means that we have a situation similar to that encountered in the one-pion 
exchange, and we could have Yukawa-like oscillations. However, some remarks are in order. At zero 
temperature and pj ~ 300 MeV we found \q \ « 200 MeV, which is the value of Aqcd- This seems to 
invalidate our one-loop calculation, since we expect non-perturbative effects to be very large under these 
circumstances. Nevertheless, pj introduces a new energy-scale into the problem (apart from Aqcd) and 
therefore the situation may be different for large values of pf. In order to investigate this possibility, we 
have calculated the evolution of go for increasing densities. The result can be seen in Fig. 5, where we 
have represented both the real part (solid line) and the imaginary part (dashed line) of qo as a function of 



Pf (all units are in GeV). As one can see, for large values of pf they become larger than Aqcd, making 
the perturbative expansion hopefully more reliable. 

A similar situation appears if one keeps the density constant and increases the temperature. In this case, 
one has to recalculate Eq. (|13| ) to take into account the temperature effects. The appropriate formulae 
can be found in ]Is[], JTsj ] , |2(fl. We address the reader to these references and will not repeat here these 
formulae. 

We have investigated the evolution of the Yukawa pole with temperature, for a fixed quark particle (minus 
antiparticle) density which corresponds to pf = 300 MeV at T = 0. This is plotted in Fig. 6, where we 
followed similar conventions to the above figure, but now in the horizontal axis one has the temperature 
(in GeV). As we found before, large values of T makes the pole appear at large momenta. We made similar 
calculations starting from larger values of pf and found no qualitative differences with the previous case. 
One expects the perturbative expansion to work better under these conditions. In fact, the asymptotic 
freedom property of QCD should manifest in this regime. 

Let us now concentrate on the physical consequences of the analytic structure of the gluon propagator 
at T = for the large-distances potential V(r). By applying the residue theorem, the pole (Yukawa) 
contribution can be written as : 

rT , . 24n(N 2 - 1) expf-^r) r , . , . . xl 
rW = ~ N (a 2 + b 2 ) r i aco < ( lrr) + bsm(q r r)} 



where 

/'(go) _ . 
= a + ib 



has to be calculated at the position of the pole . One has also to include the contribution to V(r) coming 
from the branch cut running up the imaginary axis, which we will write as Vi m (r). At zero temperature 
and large distances, this contribution can be easily approximated (proceeding as in |1C 



(UN - 2N f )g 4 {N 2 - 1) _ tt (117V - 2N f ) (N 2 - l) 
V im {r - oo) ~ - — Jjj-4-T 



with the electric mass defined as 

ml l =D{q^G) = 9 -^- 



The contribution from the branch cuts at q = ±2p/ (Friedel component) is similar to Eq. ( |ll[ ) (see 
|l0||). In Fig. 7 we have plotted the different components of V(r) (in MeV) mentioned above, at zero 
temperature and Pj — 300 MeV, as a function of the distance (in fm). The solid line gives the result of 
the numerical integration of Eq. ( [l6| ) : we will call this the exact potential. The dotted line shows the 
Yukawa component, as given by Eq. (|lj) . Also shown is the contribution from Eq. (j2Tj) (long dashes) 
and the Friedel component, as given in Rcf. [[l0| (short dashes). Finally, the dotted-dashed line is obtained 
by the sum of the Friedel, Yukawa and Eq. (p0|). As can be seen from this figure, the Yukawa component 
alone gives a reasonable approximation to V(r) for this distance range. Also, it is clear that Eq. ( |2C| ) 
overestimates the imaginary branch cut contribution for moderate distances. The Friedel component here 
is much smaller than the others. 



The situation changes for larger distances. We have represented in Fig. 8 the same pieces of the potential, 
but for a distance range r : 5 — 10 fm. Now the contribution from Vi m (r) makes the approximate potential 
closer to the exact curve. This becomes apparent when we consider even a larger r , as seen in Fig. 9. 
For this range of distances, the Yukawa component tends to disappear, while Friedel oscillations are still 
important, so that the potential can be approximated by the Friedel and Vi m (r) contributions. Anyway, 
for distances r > 5 fm, V(r) becomes too small. For the distance range where the potential is of the order 
of a few MeV, it is reasonably approximated by the Yukawa component alone. 
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IV. CONCLUSIONS 



In this work, we have considered plasma screening effects on the potential interaction for two static 'charges' in 
different models of dense matter. First, we studied the one-pion exchange in nuclear matter, as representative of 
QHD (Quantum Hadrodynamics). In addition to Friedel oscillations, which arise due to logarithmic branch cuts 
on the complex q-plane of the pion propagator, we find that this propagator shows a pole on this plane at zero 
temperature and densities above saturation. This pole gives rise to additional oscillations on the nucleon-nucleon 
potential, which we refer to as 'Yukawa oscillations'. In contrast to Friedel oscillations, which are damped as some 
power-law of the distance, Yukawa oscillations are exponentially damped with distance. Another important difference 
is the temperature behavior. While Friedel oscillations are strongly suppressed for temperatures of a few MeV, Yukawa 
oscillations are present up to T ~ 60 MeV when the density is equal to the saturation density. There is a critical 
temperature T c for each density, below which Yukawa oscillations are present. 

In order to investigate wether this phenomena are a peculiarity of this model or not, we have made further research. 
One feature of the pion exchange is the large value of the coupling constant. We considered the one-photon propagator 
in a QED plasma with an artificially increased value of the electromagnetic coupling constant. However, we found no 
Yukawa poles even if the coupling constant is taken to be of the order of unity or more. Thus, a large value of the 
coupling constant alone can not explain the origin of such pole. 

Another important feature of QHD is that it shows no asymptotic freedom. This feature originates some pathologies 
in the meson sector, such as the appearance of tachyonich branches |m| , p5| |. Therefore, one might suspect that 
Yukawa poles are also a pathological behavior of non asymptotically- free field- nuclear models. In order to investigate 
this aspect, we have studied the singlet quark- antiquark potential in a QCD quark-gluon plasma. In the case we 
studied, Yukawa poles appear, at zero temperature, when the quark Fermi momentum reaches pf > 120 MeV and 
remain for larger densities. When the temperature is increased they survive, at least up to T = 2 GeV, which is the 
maximum temperature we reached in our calculations. This situation is in contrast with the one-pion model, where 
the pole disappears at some critical temperature. It is important to notice that, although the values of q at the pole 
are close to Aqcd f° r low densities and temperatures, they become large (as compared to Aqcd) if the density and/or 
temperature increases. Since one expects to obtain asymptotic freedom in this limit, a perturbative treatment like the 
one we have considered has more chances to be valid in these circumstances. In other words, Yukawa poles in a nuclear 
medium do not seem to be associated to non-asymptotically free theories and might correspond to some underlying 
physical phenomenon. As in the case of Friedel oscillations, this could manifest in the appearance of new phases in 
dense matter. Of course, they could also be a consequence of the approximations made in our models (namely the 
one-loop approximation) and might be not present in more realistic calculations. This will be the subject of a future 
work. 

Acknowledgments 

We are grateful to Drs. V. Vento and A. Nieto for fruitful discussions. This work has been partially supported by 
Spanish DGICYT Grant PB94-0973 and CICYT AEN96-1718. 



[1] J.D. Bjorken, S. Drell, " Relativistic Quantum Fields". (McGraw-Hill, New York, 1965). S. Weinberg, "The Quantum 

Theory of Fields". (Cambridge Univ. Press, 1995). 
[2] J. Friedel, Phyl. Mag. 43 (1952), 153. Nuovo Cim.7 (1958), Suppl.2 287. 
[3] H. Sivak, Physica A129 (1985) 408. 

[4] A. Delsante and N. E. Frankel, Annals of Phys. (N.Y.) 125 (1980) 135. 
[5] T.J. Rowland, Phys. Rev. V119 N3 (1960) 900. 

[6] B.D. Serot, J.D. Walecka, "The Relativistic Nuclear Many-Body Problem", Adv. in Nucl. Phys., 16, J.W. Negele and E. 

Vogt Edts. ( Plenum, New York, 1986). 
[7] R. Machleidt. "The Meson Theory of Nuclear Forces and Nuclear Structure." Adv. in Nucl. Phys. Vol.19; J.W. Negele 

and E. Vogt Edts. (Plenum, New York, 1989.) 
[8] E. Gallego, J. Diaz Alonso, A. Perez, Nucl. Phys. A578 (1994) 542. J. Diaz Alonso, E. Gallego, A. Perez Phys. Rev. Lett. 

V.73; N.19; (1994); 2536. 
[9] J. Diaz Alonso, A. Perez, H. Sivak, Nucl. Phys. A205 (1989) 695. 
[10] J. Kapusta, T. Toimela, Phys. Rev. D37 (1988) 3731. 
[11] W. Kohn, Phys. Rev. Lett. 2 (1959) 393. 
[12] SA. Chin, Ann. Phys. 108 (1977) 301. 



8 



[13] H. Hakim, Riv. Nuovo Cimento 1 N6 (1978). 
[14] J. Diaz Alonso, Annals of Phys.; V160; Nl; (1985); 1. 
[15] T. Matsui, B.D. Serot, Ann. Phys. 144, Nl, (1982), 107. 
[16] J. Diaz Alonso, A. Perez, Nucl. Phys. A526 (1991) 623. 

[17] M.P Allendes, B.D. Serot, Phys. Rev. C45 (1992) 2975. B.D. Serot, H.B. Tang, Phys. Rev. C51 (1991) 969. 

[18] J. Kapusta, "Finite-temperature field theory". Cambridge University Press, Cambridge (1989). 

[19] K. Kajantie, J. Kapusta, Annals of Phys. 160 (1985) 477. 

[20] T. Toimela, Int. Jour, of Theor. Phys. 24, N9, (1985), 901. 

[21] C. E. Price, J. R. Shepard and J. A. McNeil, Phys. Rev. C 41, (1990) 1231 

[22] C. E. Price, J. R. Shepard and J. A. McNeil, Phys. Rev. C 42, (1990) 247. 

[23] M.J. Lighthill, "Introduction to Fourier Analysis and Generalized Functions" (Cambridge Univ. Press, Cambridge 1964). 
[24] R. J. Perry Phys. Lett. B 199, (1987), 489. 

[25] K. Wehrberger, R. Wittman and B. D. Serot, Phys. Rev. C49, (1990), 2680. 







-2p, 



2p, 



FIG. 1. Analytically-continued structure of the pion propagator (upper part) and gluon propagator (lower part) at zero 
temperature. Represented schematically are the contour paths used for integration. The Yukawa pole is indicated as an 
encircled cross. 
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FIG. 2. Second derivative of the pion polarization as a function of q/m for three temperatures at saturation density : at zero 
temperature (solid curve) the Kohn singularity at q — 2pf is apparent. As temperature increases, the curve becomes smooth. 
We showed this behavior for T = 5 MeV (short dashes) and T = 20 MeV (long dashes). 



9 



T (MeV) 



FIG. 3. Evolution of the real part of the Yukawa pole, for three values of the density, as temperature increases. The curves 
correspond to p/ = 0.3 (solid curve), pf = 0.35 (dotted) and Pf = 0.4 (dashed), which are values of the nucleon Fermi 
momentum at T = 0. 
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FIG. 4. Same as before, for the imaginary part of the pole. 
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FIG. 5. The real and imaginary parts of the Yukawa pole in the gluon propagator at zero temperature, as a function of the 
quark Fermi momentum. 



FIG. 6. Evolution of the pole of the gluon propagator as temperature changes, for a density corresponding to pf = 300 MeV 
at T = 0. 
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FIG. 7. Different components of the potential V(r) in a QCD plasma, at zero temperature and pf = 300 MeV. Solid line : 
result of the numerical integration of Eq. (|l6|). The dotted line shows the Yukawa component, as given by Eq. ( |l^ ) . Also 
shown is the contribution from Eq. ( pp| ) (long dashes) and the Friedel component (short dashes). Finally, the dotted-dashed 
line is obtained by the sum of the Friedel, Yukawa and Eq. (20). 
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FIG. 8. Same as Fig. 7, for r between 5 and 10 fm. 
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